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^ Abstract 
<' 

. Exclusive photodisintegration of the deuteron in the l-4GeV range is de- 

scribed in terms of a simple covariant and gauge invariant approach using 

> . 

OO , an effective counting rule as the hard part of the d-np vertex. At a scatter- 

\ ing angle of 9cm = 90° a prescaling behavior of the differential cross section 

\ cx with n ~ 12 is obtained; going away from 90° the value of n de- 

creases slowly, in qualitative agreement with the recent data. 
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■ I. INTRODUCTION 

The exclusive process of photodisintegration of the deuteron addresses an interesting 
interplay of nuclear and particle physics. At low energies (say, below = 0.5 GeV) con- 
ventional nuclear models based upon meson exchange which fit the NN phase shifts give a 
satisfactory description of both the energy dependence and the angular distribution of the 
experimental cross section ll|-^. However, at higher energies nuclear potential models fail 
to explain the data @J^. This is not unexpected since at Ej > 1 GeV small distances of the 
order of 0.2 fm play a role. 

Alternatively it has been attempted to describe the cross section in terms of quark-gluon 
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degrees of freedom A possible signature for the emergence of quark-gluon degrees of 

freedom would be the observation of the onset of scaling of the cross section 0. Several 
examples where this happens have been found. For example, the ^(7, 7[)n reaction above 3 
GeV appears to be consistent with counting rules at all angles [^]. 

In case of the exclusive process d{'y,p)n the dimensional counting analysis of leads to a 
differential cross section of the form 

da/dt = (1) 

where s{6) is the cm energy (angle), and n denotes the number of elementary fields in initial 
and final state (i.e. n = 13 in case of the deuteron). 

Previous data from SLAG ^ up to E^^ =3 GeV a.t 6 = 90° indicated a scaling behavior 
consistent with n = 13.1 ±0.3. More recent data from Jefferson Lab up to 4GeV confirmed 
the scaling behavior with n ^ 13 at 6 = 90° and 69°; whereas at smaller angles, 6 = 36°, 52°, 
the best fit yielded n ^ 11.5 and 11.6, respectively. This constitutes a deviation from the 
scaling behavior predicted by simple counting rules. 

A more refined approach is the reduced nuclear amplitude (RNA) approach [|7| which is also 
based upon parton exchange between the two nucleons, but takes into account some finite 
mass and higher twist effects. However, if normalized at = l.OGeV, the prediction falls 
below the data at ^ = 90° for E^ > 3GeV. On the other hand the quark-gluon string (QGS) 
model proposed in ^ and based upon Regge phenomenology appears to describe the data 
only at small t values, corresponding to small angles. 
This indicates that in practice the situation is more complex. 

The aim of this paper is to study the question of the possible origin of the apparent scaling 
and scaling violation in more detail; in particular we address whether the occurrence of 
scaling is an exclusive pQGD phenomenom or whether it can arise from different mechanisms. 
The approach in the present paper, in which the basic degrees of freedom are taken to 
be hadronic, is an extension of the one in ref. |]TD[ which predicted the cross section at 



6 = 90° for Ej > 1 GeV in fair agreement with experiment, but failed to describe the angle 



dependence. Using a simple covariant parametrization of the deuteron vertex in terms of 
a hard component (suggested by an effective counting rule) and imposing gauge invariance 
the cross section for large but not infinite s shows a "preasymptotic" scaling behavior at 
6 = 90°. This resembles the counting rule prediction, however, with an angle dependent 
value of n (which decreases away from 90°). 



II. FORMALISM 



The general covariant half off-shell d-pn vertex = m?) can be expressed as 



= T,Y + i^l^r, + ^ [T,^'^ + ^^^^U , (2) 
' 2m 2m ^ ^ 2m ^' ^ ^ 

where the Fj are scalar functions: invariant form factors. For large virtualities the contri- 
bution of the first term in eq.(0), proportional to Fi, dominates over the last three. 
Assuming that at large s only tree-type diagrams survive, the 'y+d n+p (Born) amplitude 
can be written as the sum of the pole diagrams in the s, t and u channels (Fi = G): 

+ g(-fcD7" "'""^^^y^^V ,%]C^^(nO. (3) 
s m^ 

Here s-,t-, u- variables for the d-pn vertex are: kt = [kg — q)/2, k^ = [kg + q)/2, kg = 
{p' — n')/2, p' = p + q, d' = d + q; u{p) is the nucleon spinor, E,v{d) the polarization vector 
of the deuteron {^^{d)dy = 0), C is the charge conjugation operator, (z = n,p) denotes 
the electromagnetic coupling to the nucleon, = (^2 + ^^)7'^, and F^^,^ the corresponding 
one for the deuteron (/^ = 2//^ — 1 -|- Qd)'- 
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Both 'jNN and •ydd vertices satisfy the identities: q^FJ^ = e^^, q'^Fj^^f^ = — (s — m\)gap. 
Note, s-channel accounts for the pole part of the T-matrix of the final state n-p interaction. 
In the presence of a momentum dependent vertex function the pole diagrams themselves are 



not gauge invariant. Indeed, using the Ward-Takahashi identity for the 3-point electromag- 
netic vertices, eq.@, and the Dirac equation one finds that the contraction of with the 
sum of the s-, t-, u-pole amphtudes does not vanish (i.e. the corresponding Born current is 
not conserved), if the strong d-pn vertices contain momentum dependent form factors 

Therefore the Born current in eq. @ is not complete and a contact contribution (which should 
not contain any pole- type singularities!) must be added to provide current conservation on 
the tree-level. To this end we use minimal insertion of the gauge field directly into the d-pn 



vertex [11], which gives rise to a contact amplitude, T^: 



= UdMp'h" / -Y7^{^pG{-iks - Xq/2f) + e„G(-(fc. + Xq/2)')}Cu^in'). (4) 

Jo A aQfj, 

Calculating the contraction g^T^'^ with Cp = 1 and e„ = 0, one can easily check that the total 
current is conserved |]T0|: QfiiTpoie + '^c) — 0' irrespective of the explicit form of the strong 



form factor G^—k"^), and hence the total amplitude is gauge invariant. 
To proceed we need make a choice for the vertex function G{—kf). Previous work |^ showed 
that in the energy region above > IGeV a conventional potential model wave function 
cannot describe the data. 

In the present approach it is hypothesized that the d-np vertex can be separated into two 
parts: a soft part corresponding to conventional meson exchange theory, and a small hard 
component caused by short-range phenomena. It is assumed that the soft part describes all 
low-energy (static) properties of the NN system and provides the dominant contribution to 
the normalization of the bound state wave function, while the hard part dominates the cross 
section at large virtualities. Since the microscopic structure of the short-range dynamics is 



poorly known we will use an effective counting rule prescription |T3[ to describe the hard 
part of the d-pn vertex: 



where p is the momentum of the off-shell nucleon, C is a normalization parameter and A 
is related to the inverse of the range. For the special case g = 3 [|T^ the three-pole vertex 
represents one meson propagator and two (monopole) nucleon-meson form factors. 
At the large virtualities involved obviously relativistic effects play an important role. In 
practice there exist various relativistic formulations, such as the instant form (if) formalism 
and the light-front (If) approach. Whereas in an exact calculation these are expected to 
yield the same result, this is not true in a truncated Fock space scheme. Indeed it has been 
noted ||l4| , |l6i that in lowest order (lA) the If and if approaches lead to different results. In 
particular in the lowest order Fock states in the If approach negative energy states do not 
enter. To illustrate this model dependence in the following we distinguish the covariant 
(instant form) and the light-front approach. 
Covariant Approach: Instant form kinematics 

In terms of the variables ki the vertex in eq.(^ takes on the form {i = s, t, u) 

where k'j = —al?, = —(2 — a)P, k'^ = — P, with P = s/4 — and a = 2pq/{dq) = 
1 — cos6'. Furthermore 5^ = A^/4 -|- with Oq = — my 4. 

Substituting eq.(||) in eq.(|D the cross section can be obtained straightforwardly. The abso- 
lute value of the cross section cannot be predicted and therefore the (only) parameter C is 
fixed by fitting the data at E^=l GeV. We have chosen this energy for the normalization 
since it is in this region that the microscopic d-np vertex gives a reasonable description of 
the absolute value of the empirical cross section. 

The resulting cross section s^^^ ai 9 = 90° is compared with data in fig. 1 for (7 = 3 and 
gf = 4. For comparison the results for the reduced amplitude approach of ref which is very 
close to that g = 4, are also shown. It is seen that the observed overall energy dependence 
of the cross section in the energy region 1 < < 4GeV is described well with g = 3. Only 
at the highest energies a larger value of g would fit better. 

Turning to the angular distribution one sees from fig. 2 (dot-dashed curve) that in the if 



approach the predicted angular dependence on 9 increases rapidly away from 90°. Although 
the data are available only for a few angles this is clearly not observed experimentally. 
Light-front approach 

Here the light-cone variables a and k±_ have a direct physical interpretation as the longi- 
tudinal and transverse momentum fraction carried by the nucleon in the deuteron, with 
k"^ = (s — 4m^)/4 = Q(2-"*a)^^^ • oriente the normal vector of the light-cone hyper- 

surface along q to suppress Z-graphs 0]. In general strong If form factors are functions 
of two variables, for which we can use any convenient pair from the set (a, k±, k^, k'^). For 
simplicity we will assume a factorized form of the If form factor: 

G'\a, P) = ^ with = a-^ (7) 

Here the functional dependence on k"^ is the same as in the instant form, eq.(^). The function 
0(a) in eq.(|^) goes to unity in the nonrelativistic limit. The simplest choice is k = 0, (as 
in ref |16|), but a more realistic choice is k = | [|1^ (basically corresponding to the Wick 
Cutkovsky model). At = 90° (where a = 1) both vertices are identical and the If and if 
formalisms lead to the same cross cross section (apart from a slight difference coming from 
the contact term, eq.(^). However, at 7^ 90° the form factors (6) and (7) reproduce a 
completely different dynamics. Note, that the angular dependence of the cross section arises 
mainly from the dependence of the arguments of G in eq.(^ on a = (1 — cos^). 
In G'^ in eq.(|^) these arguments are k'f = —ak"^ and k"^ = —{2 — a)k'^ . This gives rise to a 
strong increase of the cross section at forward and backward angles, which has a dynamic 
origine. Namely because of the requirement of covariance in the instant form the vertex 
dependence on a and k"^ effectively reduces to the dependence on one covariant variable 
only, say k'f (or k"^); the latter is a function of both s and 6, and hence in this case the angle 
dependence is essentially dictated by the value of the vertex parameter g in eq.(^), which 
also determines the s- dependence. 

On the other hand in the If approach the angular distribution is flatter since in this case 
the vertex depends on two variables, a and fc^, which are in principle independent. Thus, 
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a steep s-dependence of the cross section may be consistent with a flat ^-dependence. One 
sees from the curves in fig. 2 that with k = | the observed angular dependence at i^^ = 3.2 
GeV is described well. 

In fig. 3 the calculated energy dependence of the cross sections at 4 different angles is com- 
pared with the data. While the overall agreement is reasonable the data at the forward 
angles suggest a steeper increase of the cross section with energy than predicted, and more- 
over there is a discrepancy for the highest energy at 69°. Also shown are the results of the 
RNA 0] and (except at 69°) the QGS approach f^; the latter is expected to be applicable 
only at small angles. 



III. DISCUSSION 

It is of interest to discuss the underlying mechanism for scaling in the present approach. 
Using eqs.(0-§) we can express the cross section for large (but finite) s 

dt yi(s-m2)2^;-i(s-m2)3/2^^'''^^- 

For 9 = 90° one has a = 1 and f{0,s) = 1. In the relevant region of s the contribution 
of in eq.(|D is still non-negligible; for g = 3 the calculated "prescaling" behavior in the 
relevant energy region can be approximately written as ^ ~ s"""*"^, with n — 2 = 10. 
For other angles f{6, s) depends upon s through a{d, s) and therefore is a different function 
of s at different angles, effectively giving an additional power of s in eq.(8) when 6 ^ 90°. 
Therefore only for s —>■ oo the longitudinal fractions a and a' = 2 — a do not depend upon 
s: a{s — > oo) 1 — cos(6'), and one expects an (angular independent) scaling. We find that 
in the present model in the energy region 1-4 GeV n decreases slowly from n — 2 ^ 10 at 
90° to n — 2 ^ 8 at 10° (or 170°). In practice from eqs.(E|-H|) this behavior can be expressed 



as an angular dependence of the cross section on the If variable a, namely ^ oc a"^ {a > 1) 
or oc (1 — a)"^ {a < 1). 

In the past Brodsky 0] has discussed the concept of a "hadron helicity conservation law". 
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In case of the 7 + ^ j9 + n reaction it states that only hehcity amphtudes satisfying the 
condition Ap + A„ = Xd contribute at s ^ m^. Taking into account that only the Dirac 
coupling 7^) in the 7A^A^-vertex conserves hadron helicity, while Pauli coupling a^^) 
does not, we see that in the present approach asymptotic "hadron helicity conservation" will 
occur only in case the Pauli couplings are negligible in the limit s ^ m^, i.e. if the latter 
fall off at least as Note that the gauge constraint for the 3-point EM Green function, 
which does not allow any form factors in the Dirac coupling (in case of a reducible vertex 
with real photon), does not lead to any restrictions for the Pauli one. 

In practice helicity conservation predicts that for s ^ the cross section in the backward 
hemisphere receives no contribution from the neutron pole (located near 180°) but only from 
the deuteron pole which does not depend on the angle. Hence at backward angles the cross 
section would scale with s independent of angle. On the other hand at forward angles, where 
we have a competition of two different pole contributions, namely the proton one (which 
depends on 9 and s) and the deuteron one, one does not expect a unified s-dependence but 
rather an angle dependent scaling. For this reason it would be of interest to extend the 
experiment to backward angles. 

As to the differences in the results in the If and if formalisms we note that we used an 
effective counting rule for the hard part of the d-pn vertex; this is based on the assumption 
of a fixed number of constituents which is only well defined in the If approach, but not in 
the if. Therefore we consider the results of the latter (in this particular model) less reliable 
at extreme angles, which involve large t-, w-virtualities. 

We note that the deuteron vertex at large virtualities can in principle be addressed in more 
detail in semi-exclusive reactions, such as d{e, e'p)X at large in which the spectator 
proton is observed in the backward hemisphere to avoid contamination from hadronization 
products. In the past this reaction has been proposed [^] to discriminate between various 



models for the nuclear EMC effect. In the lA the deuteron hadronic tensor factorizes in terms 
of a neutron structure function and a nuclear spectral function determining the variables 
a, k± of the observed recoiling proton. Hence the cross section is directly proportional to the 
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square of the deuteron vertex. In this respect we note that the relation between the spectral 
function and the deuteron wave function is different for the If and if formalism |]15|,|T6[. (In 
fact in the covariant formulation it is not well possible to define a proper normalization.) 
This difference tends to increase with increasing a, and depend also on kj_ ||18|| . 
Finally we note that in this work we have assumed (as an extreme model) that the hard 
part of the NN interaction resides only in the (initial) deuteron vertex, and not in the final 
np state; in future work we will explore the contribution from hard contributions in the final 
np rescattering. 
Acknowledgement 

The authors are indebted to R.J. Holt, A. Radushkin and B.P.Terburg for stimulating dis- 
cussions. This work is supported in part by the Stichting voor Fundamenteel Onderzoek der 
Materie (FOM) with financial support from the Nederlandse Organisatie voor Wetenschap- 
pelijk Onderzoek (NWO). 
Figure Captions 

Fig. 1. Calculated energy dependence of ^s^^ at 9 = 90° ioi g = 3 (dashed), 4 (dotted 
curve); the solid curve is the result from Brodsky and Hiller [0]; the data are from 0,^. 
Fig. 2. Calculated angular dependence of s^^^ at = 3.2 GeV , for g = 3, in the light- 
front approach with k=0 (solid), | (dotted), 1 (dashed), and instant form (dashed-dotted 
curve). 

Fig. 3. The scaled cross section s^^^ as a function of photon energy for (a) 9 = 89°, (b) 
69°, (c) 52°, (d) 36°. The data are from ref. @J^, the solid curve shows the present results, 
the dashed curve those from (not available for 69°), and the dotted curve those from [0. 
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